Planets in triple star systems — the case of HD 188753 
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ABSTRACT 

We consider the formation of the recently discovered "hot Jupiter" planet 
orbiting the primary component of the triple star system HD188753. Although 
the current outer orbit of the triple is too tight for a Jupiter-like planet to have 
formed and migrated to its current location, the binary may have been much 
wider in the past. We assume here that the planetary system formed in an 
open star cluster, the dynamical evolution of which subsequently led to changes 
in the system's orbital parameters and binary configuration. We calculate cross 
sections for various scenarios that could have led to the multiple system currently 
observed, and conclude that component A of HD 188753 with its planet were most 
likely formed in isolation to be swapped in a triple star system by a dynamical 
encounter in an open star cluster. We estimate that within 500 pc of the Sun there 
are about 1200 planetary systems which, like Hdl88753, have orbital parameters 
unfavorable for forming planets but still having a planet, making it quite possible 
that the HD 188753 system was indeed formed by a dynamical encounter in an 
open star cluster. 

Subject headings: methods: N-body simulations - planets and satellites: individ- 
ual HD188753 - planetary systems: formation - 
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1. Introduction 

The recently discovered "hot Jupiter" orbiting the main sequence star HD 188753 A with 
period Pa = 3.3481 ±0.0009 days has significantly complicated our view of planet formation 
(Konacki 2005). The principal problem is the presence of the binary companion HD188753B, 
whose proximity severely constrains the planet formation process in the HD188753 system. 

HD 188753 was discovered and recognized as an interesting object — a stable hierarchical 
triple system — in the late 1970s (Griffin 1977). The primary HD188753A has mass itla = 
1.06 ± 0.07 M . The secondary HD188753B is itself a binary system, consisting of early 
(4.5 Gyr old) main sequence stars of total mass m B = m B1 + m B2 = 0.96 + 0.67 ± O.O7M 
(Konacki 2005). The inner (B) binary has period P B = 156.0 ± 0.1 days, semi-major axis 
a B ~ 0.67 AU, and eccentricity e B = 0.1 ±0.03 (Konacki 2005). The outer (AB) binary has 
period Pab — 25.7 yr, semi-major axis a A B = 12.1 AU, and eccentricity cab — 0.50. 

Jupiter-mass planets form in stellar and binary systems beyond the "snow line," the 
distance from the parent star where the protostellar disk is just cool enough for water ice to 
form Hayashi (1981, see also Sasselov & Lecar, 2000). For solar-mass stars, the snow line lies 
at a distance of roughly 3AU; in the particular case of interest here, it lies at r snow ~ 2.7 AU 
from HD 188753 A (Konacki 2005). Hot Jupiters subsequently sink inward due to interactions 
with the disk to a distance of ~ 10 — 20 R from the parent star. 

Truncation of the disk can prevent the formation of giant planets, and the formation 
of a gas-giant planet orbiting HD188753A in its current configuration would have been 
significantly hampered by the perturbation due to HD188753B. The disk surrounding the A 
component should have been truncated at a radius of about 1.4 AU, substantially less than 
the minimum distance (~ 2.7 AU) at which a gas giant could form Hayashi (1981) 1 . Thus 
the classical formation scenario for the planet fails, raising the question of how HD188753A's 
hot Jupiter came into being. 

Rather than exploring exotic new regimes in planet formation, this paper focuses on 
the more mundane possibility that the hot Jupiter orbiting HD188753A did in fact form 
beyond the snow line, in accordance with the standard theory of the formation of gas-giant 
planets. We propose that the present orbit came about due to a dynamical encounter in 
which the binary HD188753B became a companion of HD188753A after the planet had 
already formed and reached its current orbit close to its parent star. We present a plausible 
scenario whereby this could have occurred and argue that this is the most promising and 



1 The exact location of the snow-line and its actual relevance are discussed more extensively by Kornet, 
Rozyczka & Stepinski (2004). 
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least extreme explanation of the HD188753 system. 

The constraints on this study are simple. The planetary system (HD188753A) and 
binary (HD188753B) must have survived the interaction that brought them together, and 
the snow line in the post-encounter system must lie outside the present disk truncation 
radius, prohibiting in situ formation of the giant gas planet. We calculate cross sections for 
such encounters and discuss the corresponding formation rates in typical open star clusters. 

2. Formation Scenarios 

Given the difficulty in forming HD188753A in place, we investigate two alternative 
formation mechanisms for the HD 188753 system: 

/) The planetary system HD188753A formed as an isolated object and entered its cur- 
rent orbit around the binary HD188753B following an exchange interaction. 2 In this 
case, HD188753B must originally have had a companion star of unknown mass, which 
was ejected as a result of the dynamical encounter. There are no constraints on the 
formation of the hot- Jupiter planet, as there was no companion to perturb it. 

II) The planetary system HD 188753 A formed as a companion to another star of unknown 
mass, but in a sufficiently wide orbit to allow for the formation of the planet. An 
encounter with the binary HD188753B subsequently led to an exchange interaction, 
ejecting the unknown star and placing the binary in orbit around HD 188753 A. 

A third, but less likely, possibility is that an existing triple system was born wide enough 
to allow the formation of a planet around the primary component. Subsequent preservation 
interactions harden the outer orbit of the triple by one or more fly-bys, until the current tight 
orbit was reached. This scenario can we writen as: (AB + C — > AB + C). Our simulations 
indicate that the cross sections for this type of encounter are comparable to the other cross 
sections presented in Fig. 1 up to an orbital period of about 10 yr. The preservation cross 
section, however, continues to rise to a maximum of o ~ 7 x 10 6 AU 2 near the hard-soft 
boundary at an orbital period of about 40 yr, after which it drop sharply. This process could 
therefore have comparable weight to the other two, except that it requires a complicated 
multiple system with a planet to begin with, which ultimately may render this channel 
unimportant 



2 An exchange is one possible outcome of a three-body single-star-binary scattering, in which when the 
incoming star displaces either the primary or the secondary of the binary (see Hut & Bahcall 1983). 



-4- 



We adopt the parameters described above for the planetary system in HD 188753 A and 
the binary companion HD188753B. Both the A and B components are assumed to remain 
largely unaffected by the dynamical encounter responsible for the present system — that is, 
we assume that the interaction did not result in an encounter close enough to perturb either 
binary component. This simplifying assumption allows us to model the various interactions 
as three-body encounters, rather than considering in detail the potentially very complex 5- 
body dynamics of the entire system. We could of course relax these criteria, but we regard 
it as unlikely that such refinements will materially affect our conclusions. Since binaries and 
multiple systems appear to be the norm in open star clusters (Kouwenhoven et al 2005), 
the assumption that component A or B formed in orbit around star C does not represent a 
significant limitation on our discussion. On the basis of a limited number of 4- and 5-body 
scattering calculations, we conclude that, with the effective radii adopted below, our 3-body 
scatterings do not differ materially from their more detailed counterparts. 

We compute cross sections for the relevant exchange encounters by means of 3-body 
scattering experiments conducted using the scatter3 and sigma3 programs in the Starlab 
software environment (McMillan & Hut, 1996; Portegies Zwart et al 2001). Throughout, 
we adopt notation in which object A is HD188753A and its planet, object B is the binary 
HD188753B, and C is the unknown third star involved in the dynamical interactions respon- 
sible for the HD 188753 system we now see. The masses of the three objects are denoted 
ttia, ttib, and m c , and their (effective) radii are r^, r B , and r C - The semi-major axis and 
eccentricity of the AB binary are denoted ciab and e^s, with similar notation for the BC or 
AC binaries should they arise. 

We vary the mass mc of the third component in the initial triple between 0.1 and 
1M , and assign it an appropriate zero-age main-sequence radius. The effective radii of 
components A and B are taken to be ta = 0.05 AU and r B = 0.67 AU, ensuring minimal 
perturbation of these components during the encounters. For our purposes, a "successful" 
exchange is defined as one in which 

1. the initial conditions are stable and allow "normal" formation of a hot Jupiter, 

2. objects never approach one another within the sum of their effective radii, and 

3. the final system is a stable hierarchical system similar to HD188753, in which giant 
planet formation appears forbidden by standard planet-formation theories. 

These criteria are described in more detail in the following subsection. We note in passing 
that the choice of effective radius re is not simply a matter of our not wishing to perturb the 
B system (whose initial conditions we do not know). Our detailed scattering experiments 
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indicate that a close encounter between A or C and the tightly bound binary B can release 
sufficient energy from B to suppress the exchange process of interest here. However, we find 
that our choice of r B in the 3-body approximation is adequate to prevent this from occurring. 



In the first series of experiments we consider the possibility that HD188753A and its 
planet exchanged into an existing triple system: A + BC — > AB + C. The incomer (A) 
in our scattering calculations is HD188753A (including its planet in a 3.3 day orbit). The 
target is a triple system consisting of an inner binary star (B) and an outer component C of 
unknown mass m C - 

We vary the initial orbital semi-major axis clbc of the outer binary between 100 and 
10 6 R Q and the relative encounter velocity at infinity between lkm/s, appropriate for an 
encounter in an open cluster, and lOkm/s, to mimic a more massive parent cluster. The 
initial outer binary eccentricity esc is drawn from a thermal distribution [p(e) = 2e] be- 
tween cbc — and esc — 1) subject to additional stability constraints as described below. 
The remaining variables — the line of the ascending node, the eccentric anomaly, the orbital 
inclination angle and the moment of periastron passage — are chosen randomly, as described 
in Hut & Bahcall (1983). 

In this scenario there are no strong constraints on the initial conditions. The incoming 
star is isolated, so the presence of a planet places no restriction on its properties, and the 
target hierarchical triple system need only be dynamically stable. A coplanar prograde triple 
system (with particles 1 and 2 forming the inner orbit, and particle 3 the outer component) 
is taken to be stable if the periastron separation of the outer orbit satisfies 



(Mardling & Aarseth 1999). Here, q out = m%f (mi+m 2 ) is the mass ratio of the outer binary, 
a in is the semi-major axis of the inner binary, and e out is the outer orbital eccentricity. 
Configurations with nonzero inclination between the inner and outer orbits are expected 
to be more stable than the planar prograde case (Mardling & Aarseth 2001), so a binary 
radius based on Eq. (1) overestimates the effective size of the B binary system and hence 
underestimates the cross sections computed here. 

Some of the resulting triple systems have highly inclined outer orbits, and may be 
unstable against Kozai resonances, causing the components of the multiple system ultimately 
to merge (Ford, Kozinsky, & Rasio 2000; Lee & Peale 2003). However, we find that only a 



2.1. Insertion of the planetary system into an existing triple 




(1) 
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small fraction ( <, 10%) of our final systems could be lost by this mechanism (see also Heggie 
& Rasio 1996). 

Condition (3) above places two constraints on the orbital elements of the final system. 
First, the maximum size of a circumstellar disk after the encounter can be approximated as 

R disk ~ 0.733R> L (1 - e) L2 ° (1 + g)" ' 07 . (2) 

Pichardo, Sparke & Aguilar (2005). Here q = msl^Ai Rl is the Roche radius of component 
A, and the eccentricity dependence results from numerical scattering calculations. A binary 
born with a maximum disk radius smaller than the snow line (-Rdisk < -Rsnow) cannot form 
gas-giant planets, and our final configurations therefore must satisfy this condition. The 
second constraint is imposed by the requirement that the final triple system also be stable, 
according to Eq. 1. This criterion is least accurate for small mass ratios and inclined outer 
orbits, but we find that our cross sections are quite insensitive to the details of the stability 
criterion. 

The thick curves in Figure 1 show the results of our cross section calculations for scenario 
I, with a range of initial orbital periods and companion masses. The cross sections peak at 
initial binary periods of around 15-20 yr, largely independent of the mass of companion 
C. The cutoff at small periods stems from the stability criterion on the initial BC binary; 
systems with periods greater than ~ 10 5 years have too little energy for exchange reactions 
to form the observed AB system. 

2.2. Insertion of a close binary into a wide binary system 

In the second scenario we consider the possibility that star A and its planetary system 
formed in a wide binary system with a companion star C with orbital semi-major axis ciac 
large enough (by the criterion just described) for planet formation not to be significantly 
affected by C's presence. In this view, the binary B exchanged into this binary to form the 
observed triple system: AC + B — > AB + C. As before, we use scattering experiments to 
determine the cross section for such an exchange encounter, now subject to the constraints 
(i) that the initial AC binary has -Rdisk > -Rsnow and (ii) that the final AB binary have 
-Rdisk < -Rsnow and satisfy the triple stability criterion. For definiteness we adopt a radius of 
10 AU as being sufficiently large for a giant planet to form and subsequently migrate into a 
hot Jupiter orbit (Nelson 2000; Alibert et al. 2005). 

The cross sections for scenario II are shown as the thin curves in Figure 1. Apart from 
a slight shift to shorter periods and a reduction of a factor of roughly of two in overall 
probability, the results are qualitatively similar to those for scenario I. 
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To assess the importance of condition (3), we performed additional simulations while 
ignoring this restriction. These less restricted encounters result in cross sections which are 
larger by a factor of ~ 4 for initial periods of up to ~ 10 yr, rising to ~ 10 times larger for 
binaries with periods exceeding ~ 100 yr. The peak cross section without condition (3) is 
about (j~7x 10 6 AU 2 at an orbital period of about 40 years. 

We find that, since the vast majority of encounters are hard, subsequent interactions 
with other cluster members are unlikely to widen (or ionize) the final system again to the 
point where it no longer satisfies condition (3) — that is, where it no longer appears "unusual" 
in having a planet. Effective loss of "successful" (HD188753-like) systems is achieved only 
by coalescence during a resonant encounter, and those are relatively rare. During such an 
encounter, however, it is quite possible that all stars and the planet coalesce (Fregeau et al 
2004). 



3. Discussion and Conclusions 

We can combine the cross sections presented in Figure 1 with the binary period dis- 
tribution reported by to estimate how often encounters of the sort described here have 
occurred within ~ 100 pc of the sun, and hence how many "HD188753-like" objects one 
might expect in the solar neighborhood. The period-averaged cross section for scenario I 
is oi ~ 1.8 x 10 5 AU 2 , 3.4 x 10 5 AU 2 and 4.8 x 10 5 AU 2 for m c = 0.1 M , 0.5 and 1.0 M , 
respectively. For scenario II we find a u ^ 1.0 x 10 5 AU 2 , 2.0 x 10 5 AU 2 and 3.2 x 10 5 AU 2 
for m c = 0.1 M , 0.5 and 1.0 M . 

The overall rate for encounters of type I in a cluster of N+ stars with mean number 
density n may be expressed as 

T I = f I N*na I v. (3) 

Here // = JaJbc^ where Ja is the fraction of hot Jupiter systems like HD188753A and f B c 
is the fraction of stable triples containing a binary similar to B. We can similarly calculate 
Tj7 by substituting <7jj in 07 and /// = j 'acIb m //, where Jb the fraction of short-period 
binaries and Jac the fraction of wide binaries in which one component has a hot Jupiter. 

The values of Ja, /b, Jac an d fsc are n °t trivial to assess and therefore we retain them 
as parameters in the equations. The fraction of single cluster stars hosting a hot Jupiter 
is not well known, but we can at least estimate the other fractional parameters. The fraction 
of short period binaries is fs ~ 0.25 (Duquennoy & Mayor 1991; Pourbaix et al. 2005), 
the fraction of binaries with a hot- Jupiter in orbit around the primary is Jac ~ O.I/4, 
and the fraction of triples is fsc ~ 0.1 (Tokovinin 1999). We then estimate // ^ O.IJa, 
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fn ~ 0.025/^. 

For simplicity we assume that the secondary masses in a binary system are distributed 
with equal probability between 0.1 M and the mass of the primary, in which case the average 
cross sections become 07 ~ 3.3 x 10 5 AU 2 and 07/ ~ 2.0 x 10 5 AU 2 . The total number of 
HD188753-like systems in the solar neighborhood then is N H d ~ TN c t, where T = T/ + T n 
is the total rate for scenarios I and II, N c is the total number of open star clusters within the 
volume of interest, and T is a characteristic cluster lifetime. After substitution and scaling 
to parameters typical of open clusters, we obtain 

NHD ~ °' 23 { h [10 s AlP]") (fW?) (m) (ill) ijGyr\) " (4) 

Here n ~ 1 pc~ 3 is the stellar density of an open cluster and v ~ 1 km s" 1 is the cluster's 
velocity dispersion. Note that all clusters born over the last ~ 8 Gyr contribute, as the 
observed system is dynamically stable and the most massive star requires at least this time 
to leave the main- sequence. 

The open cluster catalog of Kharchenko et al. (2005) gives ages, distances, core radii 
and observed membership for 520 open clusters, ~ 90% of which lie within 2kpc of the 
Sun. Here we consider only the 79 clusters lying within 500 pc listed in the Kharchenko 
et al., (2005) catalog. For each of these clusters we compute the stellar density from the 
core radius and the number of stars, and the velocity dispersion (adopting a mean mass of 
0.5 M ) from the number of cluster members and the core radius from the catalog. The 
average stellar density the cluster in this sample is n ~ 6.7 stars pc~ 3 ; the mean number of 
cataloged stars per cluster is N* ~ 35. If we then adopt fjoj + fn&n — 5.3 x 10 4 AU 2 and 
sum the contribution for each of the clusters in the Kharchenko et al. catalog we arrive at 
about one planetary system with characteristics similar to HD 188753 within 500 pc of the 
Sun. 

However, the Kharchenko et al. catalog is rather incomplete. They mainly used data 
from the ASCC-2.5 full-sky survey, which is complete down to V — 11.5 mag (Kharchenko 
et al. 2005). At a distance of 100-500 pc this magnitude corresponds to that of a 0.8- 
1.5 M main-sequence star. For a Kroupa (2002) initial mass function, 87-95% of the stars 
are less massive than this, ignoring stellar remnants, and would be unaccounted for in the 
Kharchenko et al. catalog. We correct for this by increasing the number of stars in the cluster 
catalog by an order of magnitude to make up for unobserved stars. Summing the modified 
contributions for all clusters results in a total of ~ 170 planetary systems like HD188753. 

An additional correction should be made for the star clusters which are not in the 
catalog. Kharchenko et al (2005) estimate that there are about 1700 open star clusters 
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within a kpc of the sun (see also Lynga 1987; Loktin et al. 1997), whereas the Kharchenko 
et al (2005) catalog only lists 234 within that volume. Applying both corrections — for 
unobserved cluster members and incompleteness in the catalog — we estimate that, within 
500 pc of the Sun, there are about 1200 planetary systems which, like HD188753, have orbital 
parameters apparently unfavorable for forming giant planets, but which nevertheless contain 
such a planet. 

This paper was written during the MODEST-5C summer school in Amsterdam which 
was supported by the Dutch Organization for Scientific Research (NWO), the Dutch Ad- 
vanced School for Astronomy (NOVA), and the Royal Netherlands Academy of Arts and 
Sciences (KNAW). We also acknowledge the support of NASA's Astrophysics Theory Pro- 
gram under grant NNG04GL50G, and thank the other instructors at the school: Douglas 
Heggie, Piet Hut, James Lombardi, Bill Paxton, and Peter Teuben. 

Note added in proof: A recent submission of Eric Pfahl to ApJ Letters (astro- 
ph/0509490) performs similar calculations and arrives at comparable conclusions to those 
presented this paper. 
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